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^ Abstract 

Qj The cosmological origin of both dark and baryonic matter can be explained through a unified 

Q mechanism caUed hylogenesis where baryon and antibaryon number are divided between the visible 

^T") sector and a GeV-scale hidden sector, while the Universe remains net baryon symmetric. The 

"missing" antibaryons, in the form of exotic hidden states, are the dark matter. We study model- 
4^ building, cosmological, and phenomenological aspects of this scenario within the framework of 

supersymmetry, which naturally stabilizes the light hidden sector and electroweak mass scales. 
Inelastic dark matter scattering on visible matter destroys nucleons, and nuclcon decay searches 
(~| ofi"cr a novel avenue for the direct detection of the hidden antibaryonic dark matter sea. 

<N 

> 

O 

m 
m 

o 

(N 



X 



1 



I. INTRODUCTION 



The cosmological origin of both baryonic matter [T] and dark matter [2] remain an 
important mystery in our understanding of the early Universe. An array of astrophysical 
observations indicates that a fraction fib ~ 4.6% of the energy content of the Universe is 
baryonic matter, while a fraction ^dm ~ 21% is dark matter (DM) [3J. The Standard Model 
(SM) is incapable of explaining these observations, providing no viable DM candidate, nor 
a successful mechanism for generating the baryon asymmetry. Cosmology therefore requires 
new fundamental physics beyond the SM, and it is important to find ways to detect such 
new physics experimentally. 

The apparent coincidence between the densities of dark and baryonic matter, given by 
fiDivi/fib ~ 5, may be a clue that both originated through a unified mechanism. A wide 
variety of models have been proposed along these lines within the framework of asymmetric 
DM P4T6]: see Ref. [TTj for a review. In these scenarios, DM carries a conserved global 
charge, and its relic abundance is determined by its initial chemical potential. Moreover, 
if the DM charge is related to baryon number (B), then the cosmic matter coincidence is 
naturally explained for (9(5 GeV) DM mass. 

In this work, we explore model-building, cosmological, and phenomenological aspects of 
hylogenesis ("matter-genesis"), a unified mechanism for generating dark and baryonic matter 
simultaneously [El HE] . Hylogenesis requires new hidden sector states that are neutral under 
SM gauge interactions but carry non-zero B. CP-violating^ out-of-equilibrium decays in the 
early Universe generate a net B asymmetry among the SM quarks and an equal-and-opposite 
B asymmetry among the new hidden states. The Universe has zero total baryon number, 
but for appropriate interaction strengths and particle masses, the respective B charges in the 
two sectors will never equilibrate, providing an explanation for the observed asymmetry of 
(visible) baryons. The stable exotic particles carrying the compensating hidden antibaryon 
number produce the correct abundance of dark matter. Put another way, DM consists of 
the missing antibaryons. 

The minimal hylogenesis scenario, described in Refs. [T6l [18], has the following three 
ingredients: 

1. DM consists of two states, a complex scalar $ and Dirac fermion each carrying 
B = -1/2. 

2. A Dirac fermion X, carrying B = 1, that transfers B between quarks and DM through 
the gauge invariant operators [12] 

Xu%,d%dl,, X^<^ (1) 

where i,j, k label generation (color indices and spinor contractions are suppressed). 

3. An additional f/(l)' gauge symmetry that is kinetically mixed with hypercharge and 
spontaneously broken near the GeV scale, producing a massive Z'. 



^ C is charge conjugation and P is parity. 
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FIG. 1: (color online) The three steps of hylogenesis. 

With these ingredients, hylogenesis proceeds in three stages, which we illustrate schematically 
in Fig. [T] 

1. Equal (CP-symmetric) densities of X and X are created non-thermally, e.g., at the 
end of a moduli-dominated epoch when the Universe is reheated through moduli decay 
to a temperature Trh in the range of 5 MeV < Trh ^ 100 GeV ^ mx [20]. 

2. The interactions of Eq. ([T]) allow X to decay to umdR^ djik or ^'$*, and similarly for 
X. With at least two flavours of X, these decays can violate CP leading to slightly 
different partial widths for X relative to X, and equal-and-opposite asymmetries for 
visible and hidden baryons. 

3. Assuming $ and \E' are charged under U {!)', the symmetric densities of hidden particles 
annihilate away almost completely, with — )■ Z'Z' and $$* — )■ Z'Z' occurring very 
efficiently in the hidden sector, followed by Z' decaying to SM states via kinetic mixing. 
The residual antibaryonic asymmetry of $ and \E' is asymmetric DM. Likewise, the 
symmetric density of visible baryons and antibaryons annihilates efficiently into SM 
radiation. 

Both \& and $ are stable provided — m$| < nip + rrie, and they account for the observed 
DM density for + m$ ^ 5mp, implying an allowed mass range 1.7 < $ < 2.9 GeV. 

On the phenomenological side, hylogenesis models possess a unique experimental signa- 
ture: induced nucleon decay (IND), where antibaryonic DM particles scatter inelastically on 
visible baryons, destroying them and producing energetic mesons. If X couples through the 
"neutron portal" Ufidfidji, IND produces vr and f] final states, while if X couples through 
the "hyperon portal" URd^SR, IND produces K final states. These signatures mimic regular 
nucleon decay, with effective nucleon lifetimes comparable to or shorter than existing limits; 
however, present nucleon decay constraints do not apply in general due to the different final 
state kinematics of IND. Searching for IND in nucleon decay searches, such as the Super- 
Kamiokande experiment [2T] and future experiments [22] - [2l] . therefore offers a novel and 
unexplored means for discovering DM. 

Although the minimal hylogenesis model described above successfully generates the 
cosmological baryon and DM densities, two puzzles remain. Is there a natural framework to 
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consider DM as a quasi-degenerate scalar /fermion pair? Is there a mechanism to ensure the 
quantum stabihty of the GeV-scale masses for hidden sector scalars? Supersymmetry (SUSY) 
can provide answers to both questions; the DM pair ($, forms a supermuhiplet with 
B = —1/2, and the stabihty of the GeV-scale hidden sector and the ($, mass splitting is 
ensured naturally, provided SUSY breaking is suppressed in the hidden sector compared to 
the visible sector. 

The goal of this work is to embed hylogenesis in a supersymmetric framework of natural 
electroweak and hidden symmetry breaking, and to study in detail the cosmological and 
phenomenological consequences. In Section |TT| we present a minimal supersymmetric 
extension of the hylogenesis theory described above. We also address the origin of the 



nonrenormalizable nucleon portal operator X (i^^ c?^^. In Section III, we investigate 
the cosmological dynamics of supersymmetric hylogenesis, showing explicitly the range of 
masses and parameters that can explain the correct matter densities. Section IV contains a 
discussion of how such parameter values can arise in a natural way from various mechanisms 
for supersymmetry breaking. In Section |V] we investigate the phenomenology of our model, 
including IND signatures, collider probes, and DM direct detection. Our results are 



summarized in Section VI In Appendix Kl we also present an alternative supersymmetric 



model based on Higgs portal mixing. 



II. SUPERSYMMETRIC HYLOGENESIS MODEL 



In this section, we present an extension of the Minimal Supersymmetric Standard 
Model (MSSM) that can account for the dark matter and baryon densities through a unified 
mechanism of hylogenesis. In order to organize our discussion, it is useful to divide our 
model into three sectors, given by the superpotential terms 

W = VTmSSM + Wns + W^trans ■ (2) 

First, VFmssm corresponds to the superpotential of the usual MSSM with weak-scale 
superpartners; this is the visible sector. Second, we introduce a hidden sector comprised 
of new states which carry B, but are uncharged under the SM gauge group, and whose 
interactions are described by the superpotential Wns- The third term Wtrans corresponds to 
operators responsible for B transfer between the visible and hidden sectors. Baryon transfer 
operators generate equal-and-opposite B asymmetries within the two sectors, and lead to 
IND signatures in nucleon decay searches. 
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superfield 


U{1)' B R 


I'll 1 

nidden baryons Ai_2 
Yi 
Y2 


+1 —1 
-1/2 i 
+1 -1/2 i 


hidden U{1)' H 
Z' 


+1 +1 




TABLE I: New superfields in the hidden sector, with quantum numbers under B, and R- 

parity. Chiral supermultiplets Xi^2,Yi^2, H also include vector partners X^2^Yf2^H'^ with opposite 
charge assignments (not listed). 



A. Hidden sector 



The hidden sector of our model consists of (i) four vector-hke chiral superfields carrying 
nonzero 5, denoted Xi 2 and ^1^2, with charge-conjugate partners and ¥^21^ and (ii) a 
U{1)' gauge sector, with gauge boson Z' and gauge coupling e', spontaneously broken by a 
vector pair of hidden Higgs supermultiplets H,H'^. Table [T] summarizes these exotic fields. 
The superpotential is given by 



a=l,2 



+ i^xX^X: + fxyY^Y: + ixhHH'^ , (3) 



which includes Yukawa-type interactions with couplings Ci,2, Ci,2, 7, 7, and vector masses 
/ixi,2) /^yi,2 5 f^H- We also assume a canonical Kahler potential for these multiplets. Note as 
well that we have extended i?-parity to a for the Y-^'^^ multiplets. Aside from allowing 
the couplings listed above, this extension does not lead to any novel features in the present 
case, beyond those imposed by the standard -R-parity. 

After symmetry breaking in the hidden sector, the superfields Yi,2 and Yf2 mix to form 
two Dirac fermions (a = 1, 2) and four complex scalars $6 {b = 1, 2, 3, 4) with B = —1/2. 
Among these, the lightest states and $1 are stable DM. The fermionic mass terms for 
Yi^2 are (in two-component notation) 

^ferm = - Y2') My(^']+ h.C. , My ^ f "^"^^ ] . (4) 

where rj = (H), r)c ^ (H*^) are the hidden Higgs vacuum expectation values (vevs). This 
mass matrix can be diagonalized by a biunitary transformation V'^lSAyU'^ = diag(m,i,j, mijjj). 
The scalar mass terms in the basis Y = (Yi, Y2, Y^*, Y2*)'^ are 

^scalar = -YW~Y . (5) 



^ Two species Xi^2 are required for CP-violating decays (see Sec. Ill), while two species Yi,2 are needed to 
couple them to the gauge-singlet X fields. 
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The 4x4 mass matrix M~ receives contributions from F terms from Eq. ([S]), D terms, and 
soft SUSY-breaking terms 

-ifsoft 3 m^jnP + mljY^l^ + ml..\Y,'\^ + ml.\Y,'\^ (6) 
+ (feiFiFi^ + 62?2?2' + lA^Y^Y^'H + ^A^Y^Y^H' + h.c.) , (7) 



We have 



^2 ^ f Mt,My-5 + m| At \ lA,v 

also defining m~ = diag(my^, myj, = diag (my^c, myc), 5 = e'^(?7^ — 77^) x diag(0,l). 
The scalar mass matrix can be diagonalized by a unitary transformation ZM-Z^ = 
diag(m|^, ml^, m|^). Similar mass matrices arise for the X supermultiplets; for 
simplicity we assume that the fermion states Xi 2 and scalar states X12, all mass 

eigenstates. 

The U{1)' gauge sector consists of the Z' gauge boson, with mass m|, = 26'"^ {jf + r}^) , the 
Z' gaugino, and the hidden Higgsinos if, W^. The three neutralinos have the mass matrix 

M = -V2e'7] fiH ] , (9) 
V V2e'r], f^H / 

which can be brought into a diagonal form using a unitary transformation P, such that 
P^MP = diag(m^j, m^j, m^,_^). The U{iy gauge superfield mixes kinetically with the MSSM 
hypercharge, 

-^d'^J d^es'^z'^ + h.c. , (10) 

where Z'^ and Ba are the 11(1)' and U{1)y supersymmetric gauge field strengths, respectively, 
with the mixing parameter k ^ 1. 

The full particle content of the hidden sector after the spontaneous breaking of 
consists of the following mass eigenstates: three neutralinos Xii three hidden Higgs scalars 
h, H, A; two Dirac fermions four complex scalars $j; and a massive gauge boson Z'. 
The lightest Dirac fermion and complex scalar $1 are stable due to their masses and B 
charge assignments — they make up the dark matter. All other states either annihilate or 



decay into Standard Model particles as described in Sec. Ill 



Now that we have presented the ingredients for the hidden sector states, we make some 
remarks: 

• We assume that the mass scales of the hidden sector parameters lie at the GeV scale 
(with the exception of the X states). For the soft terms, this can be accomplished 



by assuming that SUSY-breaking is suppressed in the hidden sector (see Sec. IV). 
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However, the SUSY-preserving vector mass terms present a hidden //-problem; we 
ignore this issue, but in principle this can be solved by introducing an additional 
hidden singlet analogous to the NMSSM. 

• Since 2 mediates baryon transfer between the visible and hidden sectors (described 
below), IND signatures are more favorable if the DM states ($1, \E'i) are mostly aligned 
with the Yi supermultiplet. However, nonzero mixing with Y2 is induced by SUSY- 
breaking and hidden Higgs vevs resulting in a DM-Z' coupling that is essential for 
annihilation of the symmetric DM density. 

• We have imposed i? as a global symmetry. Since gravitational effects are expected 
to violate global symmetries, B violation could arise through Planck-suppressed 
operators, potentially leading to DM particle-antiparticle oscillations that can erase 
the hidden baryon asymmetry [25H28]. In our SUSY framework, these effects are 
forbidden by the extension of i?-parity. For example, the B = —1 operators 
W ~ MY^, Y1Y2H'' are allowed by U{1)' and can lead to DM oscillations, but 
they are not invariant under i?-parity. If descends from an anomaly-free gauge 
symmetry, such as U{1)b^l spontaneously broken by two units, it cannot be violated 
by gravity [29l [30] and these operators are forbidden. Thus there exists a consistent 
embedding of in a gauge symmetry that excludes the Majorana mass terms for Yi^2 
that could erase the hidden asymmetry by oscillations. 



B. Baryon transfer 



Baryon number is transferred between the hidden and visible sectors through superpo- 
tential terms Wtrans- The hidden baryon states Xi 2 are coupled to the operator U^D^Dl, 
where U^,Dj are the usual S'f/(2)i-singlet quark superfields {i,j,k label generation). We 
focus on the case involving light quarks {If^ = U^, D'^ = D^, S"^ = D^), corresponding to the 
"hyperon portal" [TU] : 

^t-ns = Yl ^^.M^aKDlS^, . (11) 
a=l,2 

with SU{?>)c indices a,/3,7 and nonrenormalizable couplings \i^2/M. Although hylogenesis 



is viable for any generational structure, Eq. (11) is the most interesting case for IND 
signatures. In contrast to non-SUSY hylogenesis models, the "neutron portal" coupling 
Xi 2U'^D'^D'^ vanishes by antisymmetry. SUSY hylogenesis therefore favors IND involving K 
final states, rather than vr, rj final states allowed in generic non-SUSY models. 

The simplest possibility to generate the nonrenormalizable coupling in Eq. ([TT| is to 
introduce a vector-like colour triplet supermultiplet P with global charges B = —2/3 and 



7 



R = 1. There are three cases to consider:^ 



A'i,2Xi,2Paf/^ + A" e^p^P^S'pD'^ + fipPaP^ (case I) 
W^trans = <( A'i,2^i,2Pa^^ + A" e^^^P^U'^S'^ + fXpP^P^ (case II) (12) 
A'i,2Xi,2Pa5^ + A" e^p^P^Df^U^ + fipP^P^ (case III) 

The SU{3)c x SU(2)l x U{1)y quantum numbers for P are (3, 1,2/3) for case I (up-type), 
and (3, 1, —1/3) for cases II and III (down-type). In all cases P'^ carries the opposite charges. 



The choice between the cases in Eq. (12) makes little difference for hylogenesis cosmology. 
However, the different cases affect the IND signals, manifested in the ratio of the rates of 
p — )■ to n — )■ channels, discussed in Sec. |v} 

Integrating out P and P'^ at the supersymmetric level generates the superpotential 



operator of Eq. ( 11 ) with Xa/M = A'^A"/ /ip together with the (higher-order) Kahler potential 



term (for case I, with similar operators for cases II and III) 

I A"P 

Including supersymmetry breaking leads to additional operators. In particular, the 
holomorphic soft scalar coupling hpPP'^ (or a squark-gaugino loop with a gaugino mass 
insertion) gives rise to the four-fermion operator Xw^d^s^ that plays a central role in IND. 



III. HYLOGENESIS COSMOLOGY 



We turn next to a study of the early Universe dynamics of our supersymmetric model of 
hylogenesis. To summarize the main ingredients: 

• We assume that the Universe is dominated at early times by a long-lived nonrelativisitic 
state <y9 (e.g., an oscillating modulus field), which decays and reheats the Universe 
before the onset of Big Bang nucleosynthesis (BBN) |20j . 

• Nonthermal CP-symmetric densities of Xi and Xi states are populated through 
decays. Depending on their specific origin, the scalar or the fermion can be created 
preferentially [3lJ. CP- violating decays of Xi and X\ generate equal-and-opposite 
asymmetries in quarks and hidden sector baryons (\l/, $), while the total baryon number 
is conserved.^ 

• A hidden f/(l)' gauge sector allows for cascade decays of heavier B = —1/2 states 
($2,\l/2, etc.) into the lightest states $ = $i and = \&i that are DM. Both states 



^ We consider the interactions of cases II and III separately, although in general both may arise 
simultaneously. The simultaneous presence of both sets of couplings leads to strangeness-violating 
interactions that may be constrained by flavor violation constraints that we do not consider here. 

^ We neglect CP- violating decays of X2 and X2, which in principle can also contribute to B asymmetries. 
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can be stable provided the condition |m,i, — m$| < mp + me is met. Also, the symmetric 
DM densities annihilate efficiently through the light Z', which decays to SM states via 
kinetic mixing with hypercharge. 

Below, we ffist compute the CP asymmetries for Xi and Xi decays. Second, we ensure 
that the successful predictions of BBN are not modified by hidden sector decays into SM 
particles. Third, we solve the system of Boltzmann equations for hylogenesis, incorporating 
all of the aforementioned ingredients, to compute the baryon asymmetries. There are 
significant differences compared to nonsupersymmetric hylogenesis [16]; in particular, the 
DM masses m$, and the ratio of $ to ^ states can be different, with implications for 
IND phenomenology. 



A. CP-violating asymmetries 



Visible and hidden B asymmetries are produced by CP violation in the partial decay 
widths of Xi 

Xi UB.dRSB. + (IrSrUr + SrUrcIr , Xi ^i^* , (14) 

due to interference between tree-level and one-loop amplitudes, shown in Fig. |2} The 
corresponding CP asymmetry is 



ex 



3 [lm(AKiC2 A2)mx, + Im(AlCi*C2A2)mx,] m^^^ 



r(Xi URdnSR) - T{Xi URd%s%) + perms. (15a) 

(15b) 



647r3M2(m2,^-m|J(|CiP + |CiP) 

We assume that Tx^ is dominated by the two-body decay to final states.^ For ex > 0, 
a positive net B asymmetry is generated in the visible sector. By CPT invariance, the 
decay rates for Xi and Xi are equal, and so an equal-and-opposite (negative) B asymmetry 
is generated in the hidden sector. Additional contributions to ex from Xi — > URdRSR, 
arising through SUSY-breaking or at one-loop, may be subleading provided squark decays 
are kinematically available. 

The baryon asymmetry can also be generated through the decays of the scalar component 
of the Xi superfield, Xi, via interference of supersymmetrizations of the diagrams in Fig. [2] 



In the supersymmetric limit, the CP-asymmetry due to Xi is equal to Eq. (15). However 



Xi decay can populate preferentially or $, due to the different hidden sector decay rates 

r(Xi ^ $*$*) = i|£mx, , TiX,-^^,^,)=^-^mx,. (16) 
167r Idtt 



In what follows, flavor indices i, j for hidden sector states are implicitly summed over in final states. 
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'■R 



X 



FIG. 2: Representative diagrams contributing to Xi — t- qiqRjQRk decays which are responsible for 
the generation of the baryon asymmetry. 

For Xi decays, the primordial ratio 

r = n^,/nq, (17) 

of charge densities $ is equal to unity. However, Xi decays can deviate from r = 1 for 
ICil 7^ ICil- As we discuss below, IND signals can be significantly enhanced if the heavier 
state is overpopulated compared to the lighter state {e.g. r ^ 1 for > m$). 

The dark matter abundance is given by 



Qb mp{l + r) 

where we have neglected the contributions from the DM anti-particles. This is appropriate 
in the limit of completely asymmetric DM populations. The allowed DM mass window, 
including the uncertainty in ^Divi/^b ~ 5, is then 

l.SGeV < mvt, m$ < 3.4GeV. (19) 

More specifically. Fig. |3] shows the allowed mass range for (blue) and m$ (red), for a 
given value of r. In the r — > {oo) limit, only $ [^>) is populated and its mass is required 
to be approximately 5mp/2 to explain the DM density; the underpopulated ^ ($) state is 



constrained within the range of Eq. (19) by the stability condition — m$| < nip + mg. 



B. Decays and Annihilations of SUSY States 

Aside from the stable DM states $i and ^^i, the hidden sector contains numerous states 
that decay, producing additional SM radiation. These decays, listed below, must occur with 
a lifetime shorter than about one second to avoid conflicts with BBN predictions. 

• The Z' gauge boson decays to SM states via kinetic mixing with the photon, requiring 



10 




Dark Matter Mass [GeV] 



FIG. 3: (color online) Allowed masses for the scalar $ and fermion ^ components of dark matter. 
For a fixed value of riij, /n$ , the shaded region shows the entire mass range of ^ (blue) and $ (red) 
that reproduces J^DM/^b ~ 5 and satisfies the stability requirement |m,i, — m$| < nie+mp. Shifting 
^Du/^h by +(— )6% moves the allowed region right (left), as indicated by the dashed contours. 

K > 10^^^ (m^//GeV)^^ ^32j, while k < 10^'^ is consistent with existing limits for 
mz' ~ GeV [MlEl].'^ 

• For the hidden Higgs states, the heavy CP-even state decays H Z'Z', while the 
CP-odd state decays A — )■ Z'h, through U{iy gauge interactions. Since the lighter 
CP-even state h is necessarily lighter than the Z', it must decay to Standard Model 
fermions either via loop-suppressed processes [35] or via D-term mixing with the MSSM 
Higgs [36] . For h masses above the two-muon threshold, the mixing process dominates 
requiring k > 10~^. 

• The heavy dark states $i [i > 1) and \&2 cascade down to $i and by emitting Z' 
and h bosons. 

• The hidden neutralinos can decay Xi ~^ $i^^i,$*\E'i provided this channel is open 
(assumed below). If this channel is closed, then the lightest state xi is stable, providing 
an additional DM component, and must annihilate efficiently via the t-channel process 
XiXi ^ Z'Z'. 



^ We assume a stronger condition k > 10 ^g^{mz' /GeV) ^{T^n/GeVY^^ sucli that tlie hidden and visible 
sectors are in kinetic equilibrium at T < Trh- |32] 
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In addition, the lightest supersymmetric particle within the MSSM decays to hidden states 
through mixing of the hidden and MSSM neutralinos induced by k. This mixing generally 
has a negligible effect on the mass eigenvalues [36]. 

The symmetric DM densities of "^i^i and $1$^ annihilate to Z' gauge bosons. In the 
case where DM is nearly aligned with the Y2 multiplet, the cross sections are given by 



e'^ / e'^ / 

{(Tv)^q,^z'Z' = iQ^^2 V ^ ~ ^I'/'^l^ {av)^^*^z'Z' = ^g^^2 V ^ ~ "^I'MI • (20) 

To have (av) ^ 3 x 10~^^ cm^/s for efficient annihilation, we require e' > 0.03 |37j . 

The presence of light hidden neutralinos allows for the chemical equilibration of baryon 
number between $ and The most important process is ^ \I'i\I'i which 

transfers the B asymmetry from the heavier DM state to the lighter state. This effect 
is phenomenologically important for IND, potentially quenching the more energetic down- 
scattering IND processes. The transfer arises from the supersymmetrization of the hidden 
gauge and Yukawa interactions, which, in the mass basis, takes the form 

if D <^i^, [(^-V2e'Z*,V,2Pik - lZ:,V,2P2k - lZ*^V,^P:,^^ Pl 

+ {V2e'z:^u,2Pl, - i*z:,u,,p;, - rz:,u;2P;,) Pr] xk + h.c. (21) 

In the limit where the dark matter is mostly aligned with the Y2 supermultiplet, the 
interaction simplifies to 

^ D v^e'$i^i {aPL + bPn) xi + h.c, (22) 

where a = —^2e'Z\^, b = \/2e'Z^2^ and Xi is the f/(l)' gaugino. Note that even if the 
mixing due to U{iy breaking can be neglected, the scalars Y2 and Y2* will still mix via the 
soft 6-term. For m$ > m,j, the s-wave contribution to the thermalized \l'i\l'i cross 

section for this interaction takes the form 



, \ _ 1 L ml 

^^''^*™"8vr(M2 + m|-m|)2\/^ x 

X (2m| [ (|a|^ + \b\^)ml + (|ap + \b\'^)iab* + a*b)m^m^ + 2\a\'^\b\'^ml ] 
-ml\{a'^ + b'^)m^ + 2abm^,\'^) . (23) 

In our numerical calculations we use a = b = e' for simplicity. In this case the transfer cross 
section reduces to 

The cross section for the reverse process \E'i\E'i — )■ is related to $i$i — )■ \E'i\l/i by the 

detailed balance condition 

{(7v),s>'f^^^ = {n^/n^f {crv)q,^^^q,, (25) 
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where the equihbrium distributions are given in Eq. (34). We discuss depletion of the 



heavier DM state in more detail below. Before doing so, however, let us mention that the 
transfer process is not generic, and may be absent in other constructions. In particular, this 
is true of the alternate Higgs portal model presented in Appendix lAl 



C. Boltzmann equations 



The generation of the visible and hidden B asymmetries during reheating is described by 
a system of Boltzmann equations: 



3Hp^ 
3Hs^ 



^ipP<p 



T 



Pv 



riB = - SHub + (exA/x + ^X'^x) 



^ippip 



(26a) 
(26b) 

(26c) 



Here, is the energy density of the modulus field ip, s is the entropy density, and ub 
is the visible B charge density (the hidden B asymmetry is —ub)- Also, Af-^ is the 
average number of Xi or its superpartner produced per modulus decay, while e^^ is the 

CP- asymmetry from Xi,Xi decay. In the supersymmetric limit ex 
decay rate determines the reheat temperature 



ejj. The modulus 



RH 



45 



.1/4 



4vr3(7,(TRH) 



r^Mpi. 



(27) 



The total modulus decay rate is given by 



47rA2 



(2^ 



where we take A = 2.43 x 10^^ GeV to be the reduced Planck constant. Along with the 
Friedmann equation = (87rG/3)(p^ + p,.), where pr = {ti"^ / 30) g^T"^, Eqs. (26) form a 
closed set and can be solved using the method of Refs. [IHl IH] • Here g^: is the energy density 
number of relativistic degrees of freedom. Instead of entropy, one can also solve for radiation 
density. We take = 2000 TeV (corresponding to Trh ~ 270 MeV), Afx = Af^ = 1 and 
ex = ^x ~ "^-^^ ^ 10~^. This decay asymmetry can be generated for example by taking 
|Aa| ~ 1, iCal ~ 0.1,^ maximal CP-violating phase and M ^ mx^ ~ mx2/3 ~ 1 TeV. 
These parameters reproduce the observed baryon asymmetry tib = ub/s ~ 8.9 x 10~^^. 
Numerical solutions to the reheating equations for these parameters are shown in Fig. |4j 
The modulus field ip decays into radiation and the heavy states X, which immediately decay 
asymmetrically in the visible and hidden sectors, generating the baryon asymmetry and the 
dark matter abundance. 



^ The magnitude of the coupling constants is chosen to be consistent with hidden sector SUSY-breaking as 



discussed in Sec. IV 
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FIG. 4: (color online) Solutions to the reheating Eqs. (26a, 26b 26c) and DM production, described 
by Eq. ([29]). Here = p^c? jm^ and Ni = rua? for i = B, ^, ^, ^>*. 



The production of dark matter and its dynamics are described by a system of four 
Boltzmann equations for $i and their antiparticles * which take the form 



Hi 



(29) 



for 2 = \E'i, $1, Here Mi^x is the average number of species i produced per X decay. 
For i = \1/ (we drop the subscript 1 from hereon) we have 



T{X ^$) + 2T{X* ^ ^^) 



X 



Similar definitions hold for \E', $ and $*, so that 

Af^+Afq>-Afq,- A/",!,. = 2ex + 2e^. 



(30) 



(31) 



The last term on the right hand side of Eq. (29) describes the production of the species i 



through modulus decay into X which promptly decays into i. The quadratic collision terms 
Ci describe the particle-antiparticle annihilations as well as the transfer reaction f-)- $$, 



We assume that the heavier dark states ^2 and <J>i, z > 1, decay to ^'i and $1 sufRciently fast. 
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FIG. 5: (color online) Solution to the Boltzmann equations for the yields Yi = m/s as a function of 
the scale factor a. The plot on the left shows the evolution for the case when the transfer reaction 
-f-T- is turned off, while the plot on the right shows the outcome when it is active. The 
transfer drives the dark matter population into lighter state, ^ in this case. The DM (anti-DM) 
abundance is indicated by solid (dashed) lines, with dark (light) lines referring to the fermion 
(scalar) component. The parameters used are described in the text. 

required by supersymmetry. The collision terms for z = ^, $ are 

Cq, = - {orv)ii,q,^z'z' [^M-^^ " {n^Y] - {crv)^^^^q, [nl - {n^/rQf n|] (32) 

C$ = - (o-t;)$$.^z'Z' - [n^Y] - (crw)$$^vtvi> [n% - {rQ/rQf n%] , (33) 

where 

n? = £-,Tm^K,{m,./T) (34) 

is the Maxwell-Boltzmann equilibrium number density for a particle of mass rrii with Qi 
internal degrees of freedom. The collision terms for the antiparticles are identical, with the 
replacements \1/ — and $—)>$*. 

The solutions to the Boltzmann equations for the yields Yi = rii/s are shown in Fig. |5] 
for = 1.9 GeV, m$ = 2.2 GeV, = 5 GeV, e' = 0.05. We consider two cases. In 
the plot on the left, we show the limit where {av),^„^,^^^J^ = 0; this can occur when the 
rate is mixing-suppressed, for a heavy gaugino, or within models with alternative symmetric 
annihilation mechanisms (see Appendix). With the transfer turned off, the scalar and fermion 
DM sectors are decoupled. The resulting DM abundances are determined by the X and X 
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Trh = 267 MeV, nou/^h = 5.0 




1 10 100 

[GeV] 



FIG. 6: (color online) The ratio rix^i/n^ for the allowed range of mass splittings Am = m$ — 
and relevant values of the hidden gaugino mass m-^ . At each point in the plane the DM abundance 
is fixed to be JloM/^^b = 5.0. Contours of constant m^r (in GeV) are also shown. The gray contour 
shows the CMB constraint for DM annihilations from Ref. |42j . Points to the right of this line are 
excluded. 

decay asymmetries. In this limit, the dark sector reduces to two independent copies of the 
standard asymmetric DM scenario. We show the case where \l/ and $ are populated equally 
by the X decays, but, in general, the asymmetries can be different for the scalar and fermion 



DM, as discussed in Sec. Ill A 



In the plot on the right we show the result when the transfer is efficient, driving the dark 
matter population into the lighter state Since the asymmetry is also transferred into the 
lighter state, the annihilation rate is enhanced, resulting in a highly asymmetric final 
abundance. The heavier state, on the other hand, freezes out with nearly equal abundances 
of particle and anti-particle, which are about an order of magnitude smaller than that of 
The transfer reaction does not affect the production of the net hidden sector baryon number 

+ n$ — — n$. . Its evolution is shown in Fig. |4j Note that 

nq, + — — n$* = 2nB (35) 

as required by B conservation in hylogenesis. 

The ratio of the abundances of \1/ to $ is important for IND. We study the effect of varying 
the mass splitting Am = m$ — rriq, and the gaugino mass on n^/n^ in Fig. [6] For each 
point in the parameter space, we solve the reheating and DM production equations and plot 
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the final value of n<^/ni^. The reheat temperature, asymmetry and gauge coupling strength 
are the same as for Fig. [5| Setting the DM abundance to the observed value fixes mx^. 
Light gauginos and large DM mass splittings make the transfer more efficient, increasing the 
abundance of the lighter state relative to the heavier one. For small Am or heavy the 
transfer rate is suppressed. 

If the symmetric density does not annihilate efficiently, residual annihilations during the 
CMB era can inject enough energy to alter the power spectrum. The WMAP7 constraint 
on the annihilation rate for Dirac fermions or complex scalars is given by [12] 

niVLifiav) 2.42x10^27^^3/ 

'^7^-^^ < T^r 36 

^^DM GeV 

where i = $ and fij/fioM is the fraction of total DM abundance in species i. This 
constraint is shown in Fig. [6] by the gray line (parameter space to the right is excluded) . 
For the parameters we have chosen, symmetric annihilation is only marginally efficient, and 
transfer processes through the light gaugino help achieve efficient annihilation (hence, the 
large region is excluded). For larger gauge coupling e' = 0.1, symmetric annihilation 
is much more efficient, and CMB constraints are evaded in the entire parameter region in 
Fig.|6} 

Finally, let us mention that we have not included any baryon washout processes in our 
Boltzmann equations. For Tjih <^ "m^Xi, ^q, the only such processes that are allowed 
kinematically are ^ 3q and the corresponding crossed diagrams. These transitions 
require the exchanges of massive intermediate P and X states to occur. These processes are 
therefore well described by effective operators of the kind 

^ ~ -^uds^<!>, (37) 

^IND 

where the scale Ajnd is defined in Section |V] (see Eq. (52)) and the order of fermion 
contractions depends on which UV completion is used in Eq. (12). We find that the 
corresponding cross sections are safely smaller than the stringent limits found in Ref. |13] 
provided mx, n^-p ^ 300 GeV. For example, we can estimate the cross section for ^$ ^ 3q 

as 

M - f-V - (4 X 10- GeV-) ry^V i^y. (38) 
^ ' \A7rJ Af^D ^ ^ V AiND ; V3 GeVy ^ ' 

The authors of Ref. [13] found that washout is negligible for 

mnMiav) < 10-'' GeV-\ (39) 



which is easily satisfied by Eq. (38). Even if these processes were important, the Xi decay 
asymmetry can be adjusted to compensate, as long as the couplings satisfy the conditions 
imposed by SUSY breaking discussed in the next section. Thus, our omission of baryon 
washout effects is justified. 
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IV. SUPERSYMMETRY BREAKING 



Our model for asymmetric antibaryonic dark matter typically requires light hidden scalars 
with masses of a few GeV to obtain an acceptable dark matter abundance. For such masses to 
be technically natural, the size of soft supersymmetry breaking felt by the light states should 
also be near the GeV scale. This is much smaller than the minimal scale of supersymmetry 
breaking felt in the MSSM sector, which must be close to or above the TeV scale to be 
consistent with current experimental bounds. 

Such a hierarchy between visible and hidden sector soft terms can arise if the hidden sector 
feels supersymmetry breaking more weakly than the visible. We examine here the necessary 
conditions for this to be the case based on the interactions required for hylogenesis. We 
also discuss a few specific mechanisms of supersymmetry breaking that can give rise to the 
required spectrum. Motivated by the desire for large moduli masses, which are frequently 
on the same order as the gravitino mass, we pay particular attention to anomaly mediation. 



A. Minimal Transmission of Supersymmetry Breaking 



The interactions we have put forward in Section [IT] will transmit supersymmetry breaking 
from the MSSM to the hidden sector. Thus, a minimal requirement for small hidden-sector 
soft terms is that these interactions do not themselves create overly large hidden soft masses. 
We begin by studying these effects. 

The states that we wish to remain light derive from the y/g^ and if'-'^-* chiral multiplets. 
These multiplets do not couple directly to the MSSM, but they are connected indirectly by 
their interactions with the X^^^ states and f/(l)' — f/(l)y gauge kinetic mixing. Thus, the 
X^^^ multiplets and the gauge kinetic mixing will act as mediators to the hidden states. 

Beginning with the X multiplets, they will feel supersymmetry breaking from their 
direct couplings to the quarks and the triplet P'-'^^ given in Eq. (12). The transmission 
of supersymmetry breaking can be seen in the renormalization group (RG) equations of the 
soft scalar squared masses of Xi'^ and P^'^) (assuming interactions as in case I): 

(47r)2— ^ = 6\X'a\'{mj,^+mlj. + ml + \Ayf)+4\Ca\\mj,^+2ml + \A^f) 



dt 

(4vr)2^ = MCaWml, + 2ml. + l^cJ') (40) 



^ ^ qo 



= 4|A'f + ml. + ml.) - ^^gl\M,\^ + . . . 

where the Ai are trilinear soft terms corresponding to superpotential operators, and we have 
dropped subleading hypercharge contributions for P^'^'^. We see that P and will typically 
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obtain full-strength (TeV) soft masses from their direct coupling to the gluon multiplet, and 
these will be passed on to X in the course of RG evolution. Recall as well that and A" 
must both be reasonably large for hylogenesis to work. 

Turning next to the Y multiplets, we find for the Yi 

a 

(41) 

a 

The Y2 multiplets are also charged under the f/(l)' hidden gauge symmetry, which mixes 
kinetically with hypercharge with strength k. This leads to additional contributions to the 
running jHHin]. At leading non-trivial order in k, we have 

-8e''(|Mf + /t'lMiH + 2e'^Sz' - 2k^ gie'Sy 

(42) 

-8e''(|Mf + ft'lMiH - 2e'^Sz' + 2k g^e'Sy , 



where Sz' = ti{Q'm'^), Sy = tr(Fm^), gi = ^/bJSgy, and Mi is the hypercharge gaugino 
(Bino) mass. The RG equations for the soft mass of H {H^) has the same form as for Y2 
{Y2) but with signs of the last two "5"' terms reversed. 

In addition to these RG contributions, there is D-term mixing between hypercharge and 
f/(l)'. After electroweak symmetry breaking in the MSSM sector, this generates an effective 
Fayet-Iliopoulos [HI HS] term in the hidden sector of the form 



V^'^{\H\' + \Y2\'- - |K,f - ^FiY , (43) 

with ^Fi = -K{gy/2e')v^cos{2l3), where v ^ 174 GeV and tan/3 is the ratio of MSSM 
Higgs vevs. This term can be absorbed by shifting the hidden-sector soft masses by mf — )■ 

The RG equations we have presented here are valid down to the scale mgoft where the 
MSSM (and X and P^^^ if their supersymmetric masses are near mgoft) states should be 
integrated out. This will generate additional threshold corrections to the hidden-sector soft 
masses. However, these lack the logarithmic enhancement of the RG contributions and are 
typically subleading. Thus, putting these pieces together we can make estimates for the 
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minimal natural values of the soft terms in the hidden sector. In terms of msoft ~ ~ TeV 
and At = In (A,,, /msoft) (where A* is the scale of supersymmetry-breaking mediation), we find 



mpi,) > msoft (44) 
"^xw ^ I A' I msoft (^^^ I (45) 




"^yW ^ ICA'Imsoft l^^^j (46) 

3/2 



m (c),mjj(c) > max <( |7CA'| msoft ( | , i^Mi ( | , \ hr^v} . (47) 



6 / \ 6 / V 2e' 

Note that \/~Ki ^ 6 for A^, = Mpi. We see that the soft masses of the Y2 and H multiplets can 
be naturally suppressed relative to the MSSM for relatively small couplings. For example, 
choosing 7 = e' = 0.05, A'l = 1, k ~ 10~^ and ( = 0.1 yields soft masses for the ¥2^'' and 
H^^'^ multiplets below a few GeV. Therefore the direct coupling of the MSSM to the hidden 
sector need not induce overly large supersymmetry breaking in the hidden sector. 



B. Mediation Mechanisms 



We consider next a few specific mechanisms to mediate supersymmetry breaking to the 
MSSM and the hidden sector that will produce a mass hierarchy between the two sectors. 
Motivated by our desire for large moduli masses, which in supergravity constructions are 
frequently related closely to the gravitino mass |m HH], the mechanism we will focus on 
primarily is anomaly mediation. However, we will also describe a second example using 
gauge mediation with mediators charged only under the SM gauge groups. 

With anomaly mediated supersymmetry breaking (AMSB) [IHl EO], the leading-order 
gaugino mass in the hidden sector is 

M' = |^m3/2, (48) 

where b' = —4 is the one-loop f/(l)' beta function coefficient. A similar expression applies to 
the MSSM gaugino soft masses, but with replaced by the corresponding factor. Based on 
this comparison, we see that a much lighter hidden gaugino will arise for small values of the 
hidden gauge coupling ^T\. For example, with MSSM gaugino masses in the range of a few 
hundred GeV, the hidden gaugino mass will be a few GeV for e' /qsm ~ 0.1, corresponding 
to e' ~ 0.05. 

The hidden-sector scalar soft masses will also be parametrically smaller than those of 
the MSSM if the corresponding Yukawa couplings are smaller as well. The explicit AMSB 
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expressions for Yi and Y2 are 



m 



m 



■3/2 



(4vr) 



a, 6 



(49) 



m: 



■3/2 



(4vr) 



f(3f -4e'^) + 2f ECa-16e'' 



We also have 



my^c, while m^c = are given by the same expressions with 7 — )■ 7 and 
(a ^ Ca- The latter point also applies to myc relative to my^ but with A'^ — )■ as well. Thus, 

we find GeV-scale soft masses for ¥2'^'' and H^'^'^ (and TeV-scale MSSM soft masses) for the 
fiducial values Ca ~ 0.1, 7 ~ e' ~ 0.05, and ^3/2 ~ 100 TeV. Note that in these expressions 
we have neglected kinetic mixing effects which are negligible for k, < 10~^ ^ ^'/gi, as we 
assume here. 



The result of Eq. (49) shows that the AMSB scalar squared masses can be positive or 
negative, depending on the relative sizes of the gauge and Yukawa couplings. This feature 
creates a severe problem in the MSSM where minimal AMSB produces tachyonic sleptons. 
We assume that one of the many proposed solutions to this problem corrects the MSSM 
soft masses without significantly altering the soft masses in the hidden sector [521 [53]. In 
contrast to the MSSM, negative scalar soft squared masses need not be a problem in the 
hidden sector due to the presence of supersymmetric mass terms for all the multiplets. In 
particular, the supersymmetric mass terms we have included in Eqs. ([3, 12) for the vector- 
like hidden multiplets can generally be chosen so that only the H and H'^ multiplets develop 
vevs. 

Let us mention, however, that supersymmetric mass terms are problematic in AMSB. In 
particular, a fundamental supersymmetric mass term Mj will give rise to a corresponding 
holomorphic bilinear soft term of size bi ~ Mim^/2- If hi ^ f^to^ti l-^iP; such a term 
will destabilize the scalar potential. To avoid this, we must assume that the supersymmetric 
mass terms we have written in Eqs. ([3|[l2| are generated in some other way, such as from 
the vev of a singlet field. ^ A full construction of such a remedy lies beyond the scope of 
the present work, but we expect that it can be achieved in analogy to the many similar 
constructions addressing the corresponding fi — Bfi problem within the MSSM [52l [53] or 
beyond \51\ . 



A second option for the mediation of supersymmetry breaking that preserves the MSSM- 
hidden mass hierarchy is gauge mediation by messengers charged only under the SM gauge 
groups im [H]. The soft masses generated in the hidden sector in this case can be deduced 
from the RG equations, up to boundary terms at the messenger scale on the order of Kmsoft, 
which are safely small. Unfortunately, the f/(l)' gaugino mass generated in this scenario 



^ Note that we could have |/ip|, ^ W3/2 without any problems. In this case, the threshold 

corrections to the light soft masses from integrating out the heavy multiplets precisely cancel their leading 
contributions from RG, leading to a zero net one- loop AMSB contribution. 
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only appears at very high loop order, and tends to be unacceptably small |15]. This can be 
resolved if there are additional gravity-mediated contributions to all the soft masses on the 
order of a few GeV. The gravitino mass in this case will be on the same order as the hidden 
states. If it is slightly lighter, it may permit the decay ^'i — )■ ip3/2+^i (for m,i,^ > m3/2+m$J. 



V. PHENOMENOLOGY 



A. Induced Nucleon Decay 



Dark matter provides a hidden reservoir of antibaryons. Although baryon transfer 
interactions are weak enough that visible baryons and hidden antibaryons are kept out 
of chemical equilibrium today, they are strong enough to give experimentally detectable 
signatures of DM-induced nucleon destruction. In these events, a DM particle scatters 
inelastically on a nucleon N = p,n, producing a DM antiparticle and mesons. For SUSY 
models, the simplest IND events are those involving a single kaon, 

^A^ ^ , ^^K . (50) 

We consider only the lightest DM states \E' = and $ = $i; the heavier states are not 
kinematically accessible provided their mass gap is larger than {niN — itlk) ~ 400 MeV. 
Both down-scattering and up-scattering can occur (defined as whether the heavier DM state 
is in the initial or final state, respectively), but up-scattering is kinematically forbidden for 
|m,5 — m$| < rriN — trk- 

Assuming the hidden states are not observed, IND events mimic standard nucleon decay 
events N — )■ Kv, with an unobserved neutrino v (or antineutrino v). Nucleon decay searches 
by the Super-Kamiokande experiment have placed strong limits on the N lifetime r for these 
modes [21] : 

t{p K^u) > 2.3 X 10^^ years , r(n K^u) > 2.6 x 10^^ years. (51) 

However, these bounds do not in general apply to IND, due to the different kinematics. For 
— )■ Ku, the K has momentum pk ~ 340 MeV. IND events are typically more energetic: 
680 ^ Pk ^ 1400 MeV for down-scattering, and px ^ 680 MeV for up-scattering (if 
allowed). The Super-Kamiokande analysis assumes: for p — )■ K~^, is emitted below 
Cerenkov threshold in water, corresponding to pk ^ 550 MeV; for p — K^, the is 
emitted with 200 < px < 500 MeV. Therefore, IND is largely unconstrained by standard 



nucleon decay searches. The limits in Eq. (51) only constrain up-scattering IND in a subset 



of parameter space, whereas down-scattering provides typically the dominant contribution 
to the total IND rate IT8l. 



For fixed DM masses, IND is eitlier bichromatic or monochromatic, depending on whether up-scattering 
is allowed or not; the range in px corresponds to the allowed mass range 1.4 < m$.* < 3.3 GeV. If other 
hidden states ^'a>2 and ^b>2 are kinematically accessible, the IND spectrum can have additional spectral 
lines. 
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FIG. 7: IND processes at leading order in chiral effective theory (left, center). Gray dot shows 
effective B transfer operator, generated by P, X exchange in our model (right). 



Next, we compute the IND rates within our supersymmetric model, starting from the 



baryon transfer superpotential in Eq. (12). The vector-like squarks P,P'^ mix through the 
SUSY-breaking term =Sfsoft ^ bpPP^ to generate the mass eigenstates Pi, 2, with masses mp^ ^. 

The leading contribution to IND arises from tree-level Pi 2 exchange, shown in Fig. ^ giving 
(in two-component notation) 



1 



off 



X 



e„^^(c?^4)«*^j)* (easel) - , , 

ea/37(4%)(^fl^i?)* (case II) , = ^ 



™° I ^ {u%d^){sl^R)<l> (case III) -=1,2 ^p.^p.^x. 



.(52) 



Here, we have neglected higher derivative terms, and Ajnd characterizes the IND mass scale. 



The different cases, corresponding to different baryon transfer interactions in Eq. (12), lead 
to different fermion contractions. 



The effective IND rate for nucleon = p, n is 

F(iV^K)=n^(HiN™' 



(53) 



where n^^^ are the local DM number densities and (o"f )ind is the IND cross section. The 
IND lifetime can be expressed as 



t{N K) 



1 



^1 + r){Vt-ou/^h)mp 



T{N^K) 2pDM[r(H 



IND 



(54) 



with local DM mass density pdm = ""^^^^ + m$n<i,, and assuming the local ratio r = n^/n$ 
is the same as over cosmological scales. The IND cross section is estimated as 



(cr^^jiND 



QCD 



10 cm-^/s X 



A 



IND 



ITeV 



(55) 



with QCD scale mqcD ~ 1 GeV.^^ For r ~ 0(1), the IND lifetime is 

/ / N _ \ -1 / \ -1 

r(Ar ^ ir) ^ lO^Vrs X ' 



IND 



10-39 cm/s 



Pdm 



0.3 GeV/cm^ 



(56) 



^-'^ The cross section (o'w)ind also depends on the DM masses m>f_$, which for the purposes of dimensional 
analysis are comparable to mgcD- 
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FIG. 8: Total IND cross section ((Tt;)iND = (cf 
with m$ = {Q,Dy[/VLTa)mp — nui/ ~ 5mp — rm^i. The effective baryon transfer mass scale is Aind = 1 



_|_ |-^^-jiV<|-->x* Q^^gj, allowed range for mii<, 



TeV. Cases I, II, III correspond to different baryon transfer models considered in Eqs. ( 12|52 ) 



which is exactly in the potential discovery range of existing nucleon decay searches, provided 
the baryon transfer scale Ajnd is set by the weak scale. 

More quantitatively, we compute (af )ind using chiral perturbation theory, which provides 
an effective theory of baryons and mesons (and DM) from the underlying quark-level 



interaction in Eq. (52), following the same methods applied to standard nucleon decay jM] 
and with additional input from lattice calculations of hadronic matrix elements [55]. We 
refer the reader to Ref. [18] for further details. Fig. [s] shows numerical results for the 
total cross section (o"f)iND, over the allowed mass range for r = 1, for the three types 
of interactions in Eq. (|52j).^^ This calculation agrees well with our previous estimate in 



Eq. (55). However, since the typical IND momentum is comparable to the chiral symmetry 
breaking scale ~ 1 GeV (i.e., where the effective theory breaks down), we regard these results 
as approximate at best. The different rates for different cases (for fixed Ajnd) satisfy 
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(57) 



as a consequence of strong isospin symmetry [18]. The kinks correspond to up-scatting 
kinematic thresholds; to the left and right, only down-scatting is allowed, while in the center 
both up- and down-scattering occur. 

More generally, we show in Fig. |9] the allowed range for p and n — )■ IND 

lifetimes. We consider masses m^^^q> consistent with ^DM/^b ~ 5, for arbitrary r in the 



The total rate n includes both and final states, and the individual channels n — ?> K'g and 



— are (approximately) half the total rate. 
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FIG. 9: (color online) Proton and neutron lifetimes for different baryon tranfer models (cases I, II, 



III) considered in Eqs. (12) and (52). Black line/gray regions show lifetime range for any r, while 



blue curves correspond to particular r values. 

range < r < oo, and we take Aind = 1 TeV. For case I (II), the allowed region is shown 
in gray, with a smaller (larger) lifetime for n — )■ than p — 7- . For case III, shown by 
the black line, the p and n IND lifetimes are equal, modulo SUSY radiative corrections and 
isospin-breaking that we neglect. Specific values for r are shown by blue curves. The solid 
blue curves show the IND lifetimes for r = 1, corresponding to the calculation in Fig. |8| For 
r 7^ 1, the IND rate can be enhanced if the heavier state is overpopulated {e.g., r > 1 for 
m$ > m$ or r < 1 for < m$); on the other hand, the IND rate can be highly suppressed 
if the heavier state is depleted and up-scattering is kinematically forbidden. The dashed 
(dotted) blue curves show the IND lifetimes for r = 5 (r = 1/5). 

We make a number of comments: 

• The IND lifetimes scale as t{N K) (x Ajj^j^. Taking Aind in the range 500 GeV — 
5 TeV corresponds to lifetimes of 10'^° — 10^^ years. Lifetimes that can be probed in 
nucleon decay searches correspond to energy scales accessible in colliders (see below). 

• Both channels p — > K'^ and n — provide complementary information, and which 
one dominates depends on the underlying heavy states mediating baryon transfer. 

• The largest IND rates in Fig. |9] correspond to ^-dominated DM (r ^ 1) with m$ > m,i,. 
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and IND is dominated by — )■ "^K down-scattering. 

Lastly, we note that while the observation of IND would be a smoking gun signal for 
hylogenesis, a nonobservation does not rule out hylogenesis as a baryogenesis mechanism. 
The IND rate can be suppressed if (i) the effective scale Aind lies beyond the TeV-scale (due 
to small couplings or large mass parameters), (ii) baryon transfer in the early Universe 
involves heavier quark flavors, and/or (iii) the heavier DM state is depleted while up- 
scattering IND is kinematically blocked. 



B. Precision Probes 

At energies well below the weak scale, the light states in the hidden sector interact with 
the SM primarily through the gauge kinetic mixing interaction. The most important effect 
of this mixing is an induced coupling of the Z' vector boson to the SM fermions / given by 

~^D-KCwQTfrZ'J , (58) 

where QJ^ is the electric charge of the fermion and cw is the cosine of the weak mixing 
angle. Direct searches for a light Z' limit kcw ^ 10^'^ for niz' < IGeV [331 EH EH] with 
significant improvements expected in the coming years [571459] . 

The dark matter states in our scenario consist of a Dirac fermion or a complex scalar 
with a direct coupling to the Z' vector. With the mixing interaction connecting the Z' to 
SM fermions, this state will efficiently mediate spin-independent elastic scattering of the DM 
states off nuclei. We estimate the cross section per target nucleon to be [I6] 

.-,..-^».„,(?f)=(^)X4)^(^)X^)'. m 

where fin is the DM-nucleon reduced mass. While this cross section is quite large, the 
masses of the DM particles in this scenario lie below the region of sensitivity of most current 
direct detection DM searches, including the specific low-threshold analyses by COGENT [60j . 
CDMS [61], XENONIO [62], and XENONIOO [63]. For a DM mass below 3 GeV, this cross 
section lies slightly below the current best limit from CRESST [61]. Proposed low-threshold 
searches for DM scattering with nuclei or electrons are expected to improve these limits [65] . 

C. High-Energy Colliders 

The new heavy states required for hylogenesis couple directly to the SM and can 
potentially be probed in high-energy colliders such as the Tevatron and the LHC. In 
particular, the effective interactions induced by the vector-like quark multiplets (P, P'^) 
can generate monojets and modify the kinematic distributions of dijets. We discuss here 
the approximate limits that existing collider data places on the masses of these multiplets, 
although we defer a detailed analysis to the future. 
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Monojet signals arising from the effective four-fermion interaction (XM^d^s^)/M^ present 
in tlie minimal hylogenesis model were considered previously in Ref. [18j. More recent 
searches for monojets by the ATLAS [66j and CMS [67J collaborations limit the corresponding 
mass scale M to lie above 0.5 — 3 TeV. Note, however, that in our supersymmetric 
formulation the corresponding four-fermion operator is only generated once supersymmetry- 
breaking effects are included. This weakens the correlation between the monojet signal and 
the operator responsible for hylogenesis, although the limit does typically force the P^^'' 
multiplets to be at least as heavy as a few hundred GeV. On the other hand, this operator 
is directly related to the IND interaction. An alternative signal that can arise directly from 
the superpotential interaction is a "monosquark" q*X final state, with the squark decaying 
to a jet and missing energy. In both cases, collider limits may be weakened through cascade 
decays in the hidden sector, which could produce additional hidden photons or Higgs bosons 
that decay to SM states. 

A second way to probe our supersymmetric UV completion of hylogenesis is through the 
kinematic distributions of dijets, which can be modified by the direct production of the triplet 
P scalars (which are i?-even). On-shell production of scalar P states via the interactions of 



Eq. (12) can produce a dijet resonance. For heavier masses, the primary effect is described 



by the non- minimal Kahler potential operator of Eq. (13), which gives rise to a four-quark 
contact operator. Studies of dijet distributions by ATLAS |68] and CMS [69l [70] put limits 
on the masses of the P^^'^ scalars of 1-10 TeV, although the specific limits depend on the 



flavour structure of the quark coupling in Eq. ( 12 ) present in the underlying theory. 



VI. CONCLUSIONS 



Through the mechanism of hylogenesis, the cosmological densities of visible and dark 
matter may share a unified origin. Out-of-equilibrium decays during a low-temperature 
reheating epoch generate the visible baryon asymmetry, and an equal antibaryon asymmetry 
among GeV-scale hidden sector states. The hidden antibaryons are weakly coupled to the 
SM and are the dark matter in the Universe. 

We have embedded hylogenesis in a supersymmetric framework. By virtue of its 
weak couplings to the SM, SUSY-breaking is sequestered from the hidden sector, thereby 
stabilizing its GeV mass scale. The DM consists of two states, a quasi-degenerate 
scalar-fermion pair of superpartners. We studied in detail one particular realization of 
supersymmetric hylogenesis, considering several aspects: 

• We constructed a minimal supersymmetric model for hylogenesis. Hidden sector 
baryons are chiral superfields X and F, with B = 1 and —1/2, respectively. The 
lightest Y states are DM, while X decays in the early Universe generate the B 
asymmetries. 

• In addition, we introduced a vector-like SU{3)c triplet to mediate B transfer between 
visible and hidden sectors, and a hidden Z' gauge boson (with kinetic mixing) to 
deplete efficiently the symmetric DM densities. 
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• We showed that hylogenesis can successfully generate the observed B asymmetry 
during reheating. We computed the CP asymmetry from X decay and solved the 
coupled Boltzmann equations describing the cosmological dynamics of hylogenesis. 

• We studied how SUSY breaking is communicated between the visible and hidden 
sectors through RG effects. We also examined predictions within an AMSB framework. 
While anomaly mediation explains the late-time reheating epoch from moduli decay, 
we have not explicitly addressed the issues of tachyonic slepton masses in the visible 
sector and the origin of SUSY mass terms in the hidden sector. 

• Antibaryonic DM annihilates visible nucleons, causing induced nucleon decay to kaon 
final states, with effective nucleon lifetime in the range 10^" — 10^^ years. DM can be 
discovered in current nucleon decay searches, and this signal remains unexplored. 

• Collider searches for monojets and dijet resonances provide the strongest direct 
constraints on our model, and these signals are correlated with IND. Lifetimes of 
Iq30 _ ]^g36 ygg^j^g correspond to energy scales Ajnd ~ 0.5 — 5 TeV that can be probed 
at the LHC. 

• DM direct detection experiments and precision searches for hidden photons constrain 
the Z' kinetic mixing, although our model remains consistent with current bounds. 

We emphasize that our specific model was constructed to illustrate general features of 
hylogenesis, and certainly there are many other model-building possibilities along these lines. 
Nevertheless, it is clear that supersymmetric hylogenesis provides a technically natural and 
viable scenario for the genesis of matter, explaining the cosmic coincidence between the dark 
and visible matter densities and predicting new experimental signatures to be explored in 
colliders and nucleon decay searches. 
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Appendix A: Alternate Higgs portal model 



We present here an alternate model for hylogenesis. In this model, hylogenesis cosmology 
is the same as described in Sec. Ill: a nonthermal abundance of Xi states arises during a low 
temperature reheating epoch and decays out of equilibrium to produce equal-and-opposite 
baryon asymmetries in the visible and hidden sectors. The main difference, compared to the 
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model in Sec. [TT| is that the symmetric DM particle-antiparticle density annihilates into light 
(GeV-scale) scalars, which in turn decay to SM fermions through mixing with the MSSM 
Higgs pseudoscalar ^4°. We do not introduce a hidden sector f/(l)' gauge symmetry, and 
this allows for a reduced field content. 

We suppose the hidden sector consists of (i) three vector-like chiral supermultiplets Xi^2 
and Y, carrying B = +1 and B = —1/2, with charge-conjugate partners X^2 ^'^i 
respectively; and (ii) a single B = chiral supermultiplet N . The Xi^2 states are coupled to 



SM quarks via Eq. (11), and are responsible for generating CP asymmetries and B transfer 



between visible and hidden sectors, as in Sec. |IT} The hidden sector superpotential is 

Whs = CaXaY^ + e NHuH^ + ^NYY"^ + f^x^^aX^ , (Al) 

with couplings Ci,2, ^, k, and where Hud are the usual MSSM Higgs supermultiplets. 
Additional terms are forbidden by introducing a global Z3 symmetry under which N and Y'^ 
carry opposite charge, while other fields are neutral. Since NHuH^ breaks this symmetry 
explicitly, we take e <C 1. The soft SUSY-breaking terms are 

-^soft 3 m%\N\^ + ml\Y\^ + m^^Yf + {A^N^ + A^NH^Ha + A^NYY'^ + h.c.) . (A2) 

We assume all couplings are real, and take < such that the scalar component of 
acquires a vev (A^) = z/, which generates mass terms for the DM states. 

The physical spectrum of the hidden sector is given as follows. From the {Y, Y'^) 
supermultiplet pair, we have a Dirac fermion \E' = {Y,Y'^'^), with mass rrix^ = |^z/|, and 
two scalar states Y, Y'^, described by the mass matrix 

- - (?• fi ( I (A3) 

V / y ^kv^ + A^v m^^ + m% j XY"^ J 

The mass eigenstates are given by $1,2, with masses 



rrii, „ 

2 



^ (2 m| + m| + ± ^ {ml - + A{ikv^ + A^vf') , (A4) 



defined such that < m\^. DM consists of the states \E' and $1. For suppressed soft 
masses (m~ ^ ^1), mixing is maximal (the relevant mixing angle is 9 ^ 7r/4) and 
< < m$2. 

The spectrum of A^ states consists of a real scalar s, pseudoscalar a, and Majorana fermion 
A^. In the limit e <^ 1, the leading 0{e) effects will be mixing between these states and the 
MSSM Higgs bosons and Higgsinos. Here, we set e = 0, obtaining 

ml = 4fcV^ + Akiy , ml = -d>AkV , m ~ = 4A;V^ , (A5) 

where the soft mass has been traded for v using the minimization condition. Vacuum 
stabilility requires ml^^ > 0. These relations can be combined: w?- = ml + |m^. 
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A light pseudoscalar (m^ < m^,, m$ J can provide an efficient channel for DM annihilation 
— )• aa and $i<l>J — )• aa, with cross sections (erf) > 3 x 10~^^cm^/s required to deplete 
the symmetric DM densities. For $i annihilation, the interaction ^ D — yci^ I'^'iP gives 

= Vl- </< - ^ X 1°-" -»'/s X (^) * (^) " . (A6) 

For \E' annihilation, t-channel ^'^^ — )• aa is p-wave suppressed; however, if s is not too heavy 
{nis < 20 GeV), it can mediate efficient s-channel annihilation: 



where the second step follows by using the mass relations given in Eq. (A5) and assuming 
vfia ~ 0. For = 3 GeV, = mg = 10 GeV, ^ = 0.1, we have {av) ^ 4 x 10"^^cm^/s. 
In both cases, as long as ^ ^0.1, the cross sections can easily be large enough to deplete the 
symmetric densities. 

Depletion of the heavier DM state, e.g. \E'^ — t- for < m^, may occur through 

t-channel exchange. We require > + mx^ ~ 5 GeV to allow for decays A^ — )• (^i^ 
(otherwise A^ would be an additional stable DM component); however A^ cannot be too 
heavy since rris < 20 GeV. For ^ — m<i>j, and neglecting terms that are 

p-wave suppressed, the leading contribution to this cross section is 

, , ^'^cos2 26' r 

^""""^ = 647rm| V ^ ~ ""^^ ' ^^^^ 

where 9 is the Y mixing angle. For = 3 GeV, = 2 GeV, nig = 10 GeV, ^ = 0.1, we 
have 

(crv) ^ 4 X 10"^^ cmVs x cos^ 29 . (A9) 

The cross section for depletion of the heavier DM component is suppressed for maximal 
mixing {9 ^ vr/4), corresponding to the case when 'm'^,m'~^ <^ m|,. 

After DM annihilation, the pseudoscalar can decay to SM fermions via mixing with the 
MSSM pseudoscalar boson through the eNHuH^ term. The decay rate is 

,2 ^,^2„ 



- - - '-y< X 1 tan^^^ 

where the A^-a mixing angle is a ~ — emgf /(2 m^o)- To avoid BBN constraints, we want 
< Is ?a 6 X 10~^^ GeV~^ For a — )■ /i/2, we have 

r(a ^ fifi)^7x 10"^ GeV x tan^ /3 sin^ a (^^^ j ■ (AH) 
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The decay rate for a — )■ ss is faster by a factor 0{3) if kinematically available. 
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